Existence Criterion of Genuine Tripartite Entanglement 



00 
O 

o 

O 

Q 

o 

m 



Oh: 

-f— > 



> 

On 
O 
O 

in 

oo 

o 



X 



Chang-shui Yu and He-shan Song* 
Department of Physics, Dalian University of Technology, 
Dalian 116024, China 
(Dated: December 30, 2008) 

In this paper, an intuitive mathematical formulation is provided to generalize the residual entan- 
glement for tripartite systems of qubits (Phys. Rev. A 61, 052306 (2000)) to the tripartite systems 
in higher dimension. The spirit lies in the tensor treatment of tripartite pure states (Phys. Rev. A 
72, 022333 (2005)). A distinct characteristic of the present generalization is that the formulation 
for higher dimensional systems is invariant under permutation of the subsystems, hence is employed 
as a criterion to test the existence of genuine tripartite entanglement. Furthermore, the formulation 
for pure states can be conveniently extended to the case of mixed states by utilizing the kronecker 
product approximate technique. As applications, we give the analytic approximation of the cri- 
terion for weakly mixed tripartite quantum states and consider the existence of genuine tripartite 
entanglement of some weakly mixed states. 

PACS numbers: 03.67.-a, 03.65.-Ta 



I. INTRODUCTION 



Entanglement is an essential ingredient in the broad 
field of quantum information theory. It is the basis of 
a lot of quantum protocols, such as quantum computa- 
tion [1], quantum cryptography [2], quantum teleporta- 
tion [3], quantum dense coding [4] and so on. It has 
been an important physical resource. Recently, many ef- 
forts have been made on the quantification of the resource 
[5,6,7,8], however, the good understanding is only limited 
in low-dimensional systems. The quantification of entan- 
glement for higher dimensional systems and multipartite 
quantum systems remains to be an open question. 

Since the remarkable concurrence was presented [5], 
it has been shown to be a useful entanglement measure 
for the systems of qubit. More interestingly, based on 
the concurrence, Valerie Coffman et al [9] introduced the 
so called residual entanglement for tripartite systems of 
qubits. The residual entanglement is independent on the 
permutation of the qubits, hence can be employed to 
measure genuine three-party entanglement, i.e. the tri- 
partite entanglement, which opens the path to studying 
multipartite entanglement. Based on the motivation of 
generalizing the definition of the residual entanglement 
to higher dimensional systems and multipartite quantum 
systems, Alexander Wong et al [10] introduced the defi- 
nition of the n-tangle for n qubits with n even, however, 
the n-tangle itself is not a measure of the n-partite en- 
tanglement. Later, hyperdeterminant in Ref. [11] has 
been shown to be an entanglement monotone and rep- 
resent the genuine multipartite entanglement. However, 
it is easy to find that the hyperdeterminant for higher 
dimensional systems and multipartite system can not be 
explicitly given conveniently. In particular, so far the 
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hyperdeterminant as an entanglement measure has not 
been able to be extended to mixed systems. Further- 
more, a new method by constructing A^-qubit entangle- 
ment monotones was introduced by Andreas Osterloh et 
al [12] for pure states to measure the n-partite entan- 
glement, however, it is only confined to the systems of 
qubits and seems to be very difficult to extend to the 
case of mixed states analogously to Ref. [11]. 

In this paper, we introduce a new approach to gen- 
eralize the residual entanglement for tripartite systems 
of qubits to the tripartite systems in higher dimension. 
One knows that the key to obtaining the explicit tabc 
in Ref. [9] is the analytic expression of the concurrence 
in mixed systems of qubits. However, so far no one has 
been able to obtain an analytic expression of concurrence 
(or concurrence vector) for higher dimensional mixed sys- 
tems, which means that the expectable results for higher 
dimensional systems seems not to be obtained from the 
similar method to that in Ref. [9]. Hence, we provide 
an intuitive mathematical formulation to generalize the 
residual entanglement according to the tensor treatment 
of tripartite pure states presented in Ref. [13]. A distinct 
characteristic of the present generalization is that the 
formulation for higher dimensional systems is invariant 
under permutation of the subsystems (i.e. the qudits), 
hence can be employed as a criterion to test existence of 
the genuine tripartite entanglement (also called tripar- 
tite entanglement for convenience in the paper). Further- 
more, the formulation for pure states can be conveniently 
extended to the case of mixed states by utilizing the kro- 
necker product approximate technique [14,15]. However, 
it should be noted that the formulation is not an en- 
tanglement measure except that for tripartite systems of 
qubits due to the variance under local unitary operations. 
As applications, we give the analytic approximation of 
the criterion for weakly mixed tripartite quantum states 
(quasi pure states) and consider the existence of tripartite 
entanglement of some quasi pure states, which shows that 
our criterion can be conveniently applied in these cases. 



The paper is organized as follows. Firstly, we give the 
intuitive generalization of the residual entanglement for 
pure states; secondly, we extend it to mixed states and 
discuss the existence of tripartite entanglement of some 
quasi pure states; the conclusions are drawn in the end. 



II. EXISTENCE CRITERION OF TRIPARTITE 
ENTANGLEMENT FOR PURE STATES 

The residual entanglement for tripartite systems of 
qubits or tabc (i- e - the tripartite entanglement mea- 
sure) is given by 



t(\ iPabc)) = Vdet R=\d 1 - 2d 2 + M 3 \ , (1) 

where a constant factor is neglected and the element R^ 
of the 2x2 matrix R is defined by 

Rij — ^ ] Q,kljQ J mni^mp^nq®pqr®8tr€sk€th (^) 

with the sum being over all the repeated indices, eoi = 
-eio = 1 and e o = -en = 1; 

7 22|22i22i22 
"1 = a 000 a lll +" a 001 a 110 + a 010 a 101 + a 100 a 01l! 



<^2 = a ooo a m a oii a ioo + aooo a m a ioi a oio 

+<2000 a lll a 110 a 001 + a 011 a 100 a 101 a 010 
+ a 011 a 100 a 110 a 001 + a 101 a 010 a 110 a 001> 

ds = aoooanoflioiaon + «iiiaooiaoioaioo- (3) 

What's more, the a terms in above equations are the 
coefficients in the standard basis defined by \ipABc) = 
J2i j fc=o a v k Ni^) abc- ^ s mentioned in Ref. [9], the ex- 
pression of t(\iPabc}) can be mentally pictured by imag- 
ining the eight coefficients attached to the corners 
of a cube. The picture yields that r is invariant un- 
der permutations of the qubits, because a permutation 
of qubits corresponds to a reflection or rotation of the 
cube. It happens that the picture is consistent to the 
tensor cube introduced in Ref. [13]. In other words, a 
tensor cube of \tpABc) corresponds to a tripartite entan- 
glement measure t(\iPabc))- For convenience, we employ 
/(IV'abc)) = \t(\iPabc))\ 2 to measure tripartite entan- 
glement, which is equivalent to t(\iPabc)) from the view- 
point of entanglement measure. Obviously, f (\iPabc)) 
has the same properties to t(\iPabc})- 

According to Ref. [13], a tripartite pure state in 
any dimension can be regarded as the tensor grid 
which includes tensor cubes. E.g. let \4>abc) — 

E!j=o ELo a vk \ijk) AB ci the tensor grid of \4>abc) can 
be pictured as figure 1, which includes three tensor cubes. 
In this sense, one can draw a conclusion that tensor cube 
is the unit of tensor grid. Since every tensor cube in a 
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FIG. 1: The tensor grid of the coefficients of a tripartite pure 
state in 2 x 2 x 3 dimension. 



tensor grid can be considered as an non-normalized tri- 
partite pure state of qubits, one can get that every unit 
corresponds to the tripartite entanglement measure of 
the non-normalized pure state. Namely, the tensor cube 
corresponds to the minimal unit of describing the tripar- 
tite entanglement. Therefore, whether there exist some 
genuine tripartite entanglement can be determined by all 
the minimal units. 

Theorem l:-For any a tripartite pure state Y$>) which 
includes M minimal units mentioned above, let the 
the non-normalized tripartite pure state of qubits cor- 
responding to the ith unit be denoted by \<fii), then the 
corresponding tripartite entanglement can be given by 
f(\<Pi))- Define 



M 



£/(l<^», 



(4) 



for the state \^>), then if there does not exist genuine 
tripartite entanglement in 1^), F^)) = 0. 

Proof. It is obvious that F(\^)) = means that 
f (\(fi)) = holds for all ipi, vice versa. Since the ten- 
sor cube corresponds to the minimal unit of describing 
the tripartite entanglement, F(\^)) = shows that there 
does not exist genuine three-party entanglement in |^). 
That is to say, F 1 (| , I')) can effectively test the existence of 
tripartite entanglement in Furthermore, a permuta- 
tion of qudits corresponds to a reflection or rotation of 
the tensor grid, which is similar to that in Ref. [9] , hence 
all the tensor cubes in the tensor grid are invariant ex- 
cept the relative positions in the grid. Namely, F 1 (| I I')) is 
invariant under permutations of the qudits. 

Considering the matrix notation of \ty) = 



i=o 2^7=o 2—i 



"3-1 

fe=0 u v k 



ijk), F(\fy}) can be ex- 
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pressed as the function of |^), i.e 
F(|*» = 



JVi N 2 N 3 
\ a=l 0=1 7=1 



where AT p = 



with p = 1,2,3; s g , q = a,/3,7, 
denotes 2 x n p matrix with p corresponding to g. If the 
generator of the group SO(n p ) is denoted by S p , s q can 
be derived from \S P \ by deleting the row where all the 
elements are zero, where | | denotes the absolute value of 
the matrix elements. 

Because eq. (2) can also be written in the standard 
basis by 



Rij = ^](aoojaii r + aujCtoor — aoijdiOr — awjdoir) 

r=0 

'l a 00i a llr ~r a lli a 00r — a 01i°10r — a 10i a 01rJi l D ) 

and /(IV'abc)) = deti2, F(\^f)) can be expanded by 

N\ N 2 N 3 
a=l /3=1 7=1 

1 

E K**l l|0» K ® (Wl Sa/J 7 

fc=0 
1 

XEK*I 5 «/37 H 1 )) K ® ((fell ^7 I**) T 

fc=0 
1 

^/3 7 l|o» K ® <0 S a /J 7 I*) 

fe=0 

x <*| K <8>Vy)((k\\S a(h |**)) 

i 

x E«**l ^7 II 1 )) K ® *„) s afh I*) 

fe=0 

x (*| S T a01 1|0» (a y ® *„) S a /J 7 |**))]} 1/4 , (7) 

1 \ MO 



where 5 a/ 3 7 = s Q <£> S/3 <g> s 7 , ||0)) = 
l|l» = 



1 
1 



J ™ \ 1 
1 ) ( 1 )' and ^ " = 



(|| ))) and the superscript T denotes transposition op- 
eration. Note that Ffltf)) = i(F(|*)) + F(|*))*). Al- 
though the expanded F(|\P)) is a bit tedious, it is impor- 
tant for the extension of F(\^f)) to mixed states. 



III. EXISTENCE CRITERION OF TRIPARTITE 
ENTANGLEMENT FOR MIXED STATES 

A. Kronecker product approximation technique 

We first introduce the kronecker product approxima- 
tion technique [14,15]. For any a matrix M — [ray], with 



entries my, defined in <8> Cd 2 , M can be defined [16] 

by 



M = V^{MV£) 



R\T 2 



(8) 



where the superscript T 2 denotes partial transposition 
on the second space [17], V^'* are left ( ri g ht ) hand side 
swap operators defined as V 12 = J2ikj>k> Sjk'Sfk \ j) ® 
\k) (k'\, j,k' = l,---,d 2 ,f,k= l,---,di. The right hand 
side swap operator is defined in Cd 1 <g> d 2 and the left 
one is defined in Cd 2 <E)Cd 1 - Furthermore, Vfe — (V^ ) T = 
(V^) -1 . If d\ = d 2 , = Vjf . M has the singular value 
decompositions: 



M 



t/SV f = 



(9) 



where Uj, vi are the ith columns of the unitary matrices U 
and V, respectively; £ is a diagonal matrix with elements 
(Tj decreasing for i = 1, • • -,r; r is the rank of M. Based 
on Ref. [14,15], M can be written by 



(10) 



i=l 



folv*, where 



with Vec(Xi) — ^falui and Vec(Yi) 

Vec(A) = [aii,---,api,ai2,---,a P 2,---,aig,---,apg] T , (11) 
for any a p x q matrix A = [a%j] with entries a%j [18]. 

B. Extension of existence criterion to mixed states 

Consider F(\$>)) of pure states, the corresponding 
quantity of mixed states p is then given as the convex 
of 

F(p)=miJ2p i F(\^ i )) (12) 

i 

of all possible decompositions into pure states |<J/i) with 



p= , £pi\*i)m,Pi>o. 



(13) 



F(p) vanishes if and only if p does not include any gen- 
uine three-party entanglement. According to the matrix 
notation [7] of equation (13), one can obtain p = SffW^ , 
where W is a diagonal matrix with Wu — pi , the columns 
of the matrix ^ correspond to the vectors |^i). Due to 
the eigenvalue decomposition: p = $M&, where M is a 
diagonal matrix whose diagonal elements are the eigen- 
values of p, and $ is a unitary matrix whose columns are 
the eigenvectors of p, one can obtain ^W 1 / 2 — ^M 1 / 2 ?/, 
where U e C rx N is a Right-unitary matrix, with N and r 
being the column number of iff and the rank of p. There- 
fore, based on the matrix notation and cq. (7), cq. (12) 
can be rewritten as 

JV 

F(p) = inf E([(^ T ® U ^ ® UT ® U ^) 



4 



x A <*Pi {U®U*®U® £T)]^) 1/4 , (14) holds for any j, with (Cj)., (Cj) . denned in C dxd , a'j the 



a/37 



where 



A aPl = \ ( P 1/2 Y S T af3j \\L)) E y S a/37 (p 1 / 2 ) 

defined in C dxd <S> C^xd ® C d xd ® C<jxd, and p is defined 
in C dxd , with 

pV2 = ^$M 1/2 ) T ®($M 1/2 y®($M 1/2 ) T ®($M 1/2 y , 



||L)) = ||0011)) + ||1100))-||0110))-||1001)) 
||i2)) = («00|| + «ll||)®(«00|| + «ll||) ) 



corresponding singular value and (C 7 -) i = y ( cr j) i (^j)j- 

Analogously, r" is the rank of Bj. Hence, eq. (16) can 
be rewritten by 



1/4 



N 



(17) 

The infimum can be employed to test the existence of 
tripartite entanglement of p. 

In terms of the Cauchy-Schwarz inequality 

1/2 / \ 1/2 



Ey 2 > E^2/i and E 2 Nil > lEi^il, 



can also be expressed by 



JV 



and 



S a /3 7 — ®j = iSa/37> 



||afecd)) = ||a))<8> ||6))(8i ||c))(gi ||d)>. 



If the former two subspaces and the latter two ones are re- 
garded as a doubled subspace, respectively. E Q/ 3 7 A aP-y 
can be considered to be defined in C d 2 Xd 2 ® C d 2 Xd 2. It 
is easy to find that Ea/3 7 ^a/3 7 is invariant under the ex- 
change of two doubled subspaces. Hence, based on the 
kronecker product approximation technique, E Q/ 3 7 ^-a/s-y 
can be written by 



a/?7 i i 



(15) 



with £?• defined in C dxd <8> Cdxd and a-j the corre- 
sponding singular value. Bi = sJoiB^ which can be ob- 
tained following the procedure in above subsection is not 
given explicitly. Furthermore, r' is the rank of the ma- 
trix Ea/3 7 A a f3<y defined in above subsection. Due to eq. 
(15), eq. (14) can be rewritten by 



N 




Ul)Bj(U®U% 



It is also obvious that A n n~, is converted into A* 



1/4 



(16) 
if 



the former two subspaces and the latter two ones are ex- 
changed simultaneously. Based on the kronecker product 
approximation technique again, one can obtain that 



F(p)>inf]T 



U 1 



(18) 
1, and 
1. 



where Zj = exp(i</5j), with a;j > 0, Ej ^ 

= y jm exp(i(p jm ), with y im > 0, J2 jm y] m 
Eq. (18) has the similar form to that in Ref. [7], 
even though it is a little more complex. Therefore the 
infimum of eq. (18) can be given by max Ai(z, Z) — 

z.Z 

E»>i Aj(z, Z), where Aj(z, Z) are the singular values of 
(Ej EL z i ' ^j'm (Cj) m ) in decreasing order [7], with 



and Z=[Zn,Zi2,- • -,Zi r »,Z; 



z = [zi,z 2 
•, Z r / r "] . 

In terms of the inequality E™=i \ x i\ > £ (E"=i 
a more analogous form about eq. (17) to that in Ref. [7] 
can be given by 



21- 



/ 1 \ 1/4 JV 

^)>infU £ 



17 



1 4 (19) 

The infimum can be obtained by (p-) ^ 4 x 
^ i>1 Ai(Z)^j, where Aj(Z) arc the 

^ (Q) m ) n 



maxAi(Z) 
z 



singular values of J2 m z jrn (Cj) m ) in decreasing 

order. Both the two cases can provide the necessary 
condition for the existence of tripartite entanglement 
of a mixed state, but the sufficiency of them may be 
different. 

What's more, compared with the procedure in Ref. [8], 
it is very possible that [(Cj-) m ] m corresponding to the 

maximal ^J(<Jj) m can give the main contribution to the 
infimum of eq. (18). That is to say the lower bound of 
F{p) can be given by Ai — Ei>i ^» with Xj the singular 
values of [(C,) 1 

J 'ml max 



■5 



C. Examples 



^ii'ii ■ Therefore, one can have the approximation 



In above subsection, we have provided three different 
lower bounds for any mixed state, which can be employed 
as necessary conditions to test the existence of tripar- 
tite entanglement in principal. However, by analysis, one 
can find that the numerical realization to calculate the 
bounds for a mixed state p requires the eigenvalue decom- 
position of a matrix defined in the same dimension to that 
of <8>f = iP , which reduces the efficiency of calculation. In 
order to avoid the similar problem, an analytic approx- 
imation method was introduced for quasi pure states in 
Rcf. [19]. By utilizing the analogous method, one will 
find that eq. (17) can be simplified significantly, hence 
our criterion can work well for quasi pure states. Before 
the examples, we firstly give the analytic approximation 
ofcq. (17). 

Let J2 a f3j A a p-y in eq. (14) be denoted by A. Analo- 
gous to Ref. [19], the tensor A can be obtained by 

Am,jk 
^■I'm'J'k' 
iVi N 2 N 3 
= E EE yUlUmUjUkUl'Um/Uj'Uk, 
a=l /3=1 7=1 
1 

x Sl Pl ||0» K ® a v ) ((i\\S a01 

x (* m | sl 01 ||o» K ® *v) ms afh \KA) 

1 

x E«^*l ^ II 1 )) K ® °y) Ms afh l*i'> 

i=0 

x^ k \Sl M \\l)){<j y ®a y ){(j,\\S aPl \n>)) 
i 

-E«**l ^7 II 1 )) K ® "„) Sa/J 7 |*I'> 



^W.j'fc' - T l m Ti, ml T 0k T*, kl With T im = 



i/m,ll 
1 11,11 



(<7 !/ (g)<Tj / )((i|| < S' a/ 3 7 |*^)) 



X E(<*;l S l?1 H°)) K ® "„) S a /J 7 |*i'> 

x (tt fc | S^ 7 ||1)) (tr„ ® <j y ) ((i|| 5 Q/37 |**,))], (20) 

where denotes the ath eigenvector and all the other 
quantities are defined similar to those in eq. (7). Ac- 
cording to the symmetry of A and the kronccker product 
approximation technique in above section, A can be for- 
mally written as 

\lm,jk \ ' rpa (rpa \* rpa (rpa \* 

A l'm'.i'k' — / J 1 lm\ 1 l'm') 1 jk{- L j'k') ■ 



The density matrix of quasi pure states has one single 
eigenvalue \i\ that is much larger than all the others, 
which induces a natural order in terms of the small eigen- 
values ^ , i > 1 . Due to the same reasons to those in Ref. 
[19], here we consider the second order elements of type 



K;n)' 



In this sense, eq. (17) and eq. (18) can be simplified 
significantly: 

F(p)cF a (p)=m^\U T rU\... 

i 

F a (p) can be given by 

F a (p) = max{Ai - A 4 ,0}, 

i>l 

where Xi is the singular value of r in decreasing order. 
The tripartite mixed states introduced in Ref .[20] 



p(x) = x \GHZ) (GHZ\ + (l-x)/2(\W) (W\- 
where 

\GHZ) = -j=(\000) + \lll)), 



|W> = — (|001> + |010> + |100» 



\v)(w 



)> 



^ = —(|110) + |011) + |101)), 

can be considered as a quasi pure state for x > 1/3. 
F a (p(x)) is shown in Fig. 2, which indicates the consis- 
tent conclusion to that in Ref. [20]. What's more, for 
the quasi pure states generated by the mixture of max- 
imally mixed state(identity matrix) and tripartite GHZ 
state (The cases in 3 x 3 x 3 dimension is included.), the 
corresponding F a (p)s can all be shown to be nonzero. 
Since F a (p) is not monotone in higher dimension, the 
corresponding figures are not given here. 



IV. CONCLUSION AND DISCUSSION 

In summary, we have introduced an intuitive mathe- 
matical formulation to generalize the original tripartite 
entanglement to higher dimensional tripartite systems 
according to the tensor treatment of a tripartite pure 
state. A distinct characteristic of the present generaliza- 
tion is that the formulation for higher dimensional sys- 
tems is invariant under permutation of the qudits. When 
the formulation is reduced to tripartite systems of qubits, 
there exists an exponent \ different from the original 
one, but the change of exponent provides convenience for 
the generalization to mixed states. The formulation for 
pure states can be conveniently extended to the case of 
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x 

FIG. 2: F a (p) of the mixed state p(x) = x \GHZ) (GHZ\ + 
(1 - x)/2(\W) (W\ + |w^> (w\) vs x, x G (1/3, 1]. 



conditions to test the existence of tripartite entangle- 
ment, but the sufficiency of them may be different. How- 
ever, because the dimension of A a py is much higher than 
that corresponding to bipartite entanglement, it seems 
to be a bit difficult to directly apply to test the exis- 
tence of tripartite entanglement of a general quantum 
mixed state. Fortunately, for the weakly mixed states, 
i.e. quasi pure states, one can find that our criterion can 
be conveniently applied and is even a sufficient condition 
for the existence of tripartite entanglement. In particu- 
lar, our criterion can provide an analytic approximation. 
Since the 3-tangle is an entanglement measure, F a (p) is 
not only an existence criterion, but also an effective tri- 
partite entanglement indicator. Even though there exist 
some questions left open, the intuitive mathematical for- 
mulation of tripartite entanglement and the convenient 
extension to mixed states will play an important role in 
the further understanding of multipartite entanglement 
measure. 



mixed states by utilizing the kronecker product approxi- 
mate technique. We have presented three different lower 
bounds for F(p) of mixed states. The forms of the three 
results for mixed states are similar to those of bipartite 
entanglement [7,8]. All of them can provide necessary 
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